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SUMMARY  OF  RESULTS 


Basic  Formulae 

The  degradation  in  the  performance  of  the  Artificial  Ionospheric  Mirror  (AIM) 
due  to  phase  fluctuations  induced  by  irregularities  in  the  electron  density  of  the 
AIM  plasma  were  studied.  The  general  theory  was  formulated  on  the  basis  of 
the  path  integral  techniques  [Feynmann  and  Hibbs,  1965;  Dashen,  1979].  It  was 
found  that  the  presence  of  irregularities  in  the  reflector  reduces  the  signal  by 

L  =  -4.3  K  <  Sjipath)  >]  dB  (1) 

where  Iq,  is  the  free  space  wavenumber  and  Si  is  the  variation  in  the  eikonal 
induced  by  the  presence  of  irregular  structure  in  the  reflection  layer.  Equation  (1) 
when  applied  to  specular  reflection  from  a  rough  surface  with  random  inhomo¬ 
geneities  having  RMS  height  Ah,  gives 

C(rough  surface]  =  —  I7(hoAh)^  cos^  ^  dB  (2) 

This  result  is  consistent  with  the  one  derived  by  other  techniques  [Beckmaim 
and  Spizzichino,  1963;  Ogilvie,  1991]  and  reproduces  the  Raleigh  criterion.  For 
the  AIM  case  and  assuming  a  linear  profile  with  gradient  length  L,  the  loss  is 
given  by 

C{AIM)^  -A.ZgdB  (3) 

For  isotropic  fluctuations  and  reflection  at  an  angle  6, 

OO 

g  =  2n{koL)(^\n  J  Sn{k)kdk  (4) 


V 


where  Sn(k)  is  the  spectrum  of  the  electron  density  fluctuations.  For  a  Gaussian 

spectrum  with  correlation  length  £o. 

9  [^ ,  ,  .  ^  .  r  ^ .  1  +  COS  0 

9  =  ^olo){i^oL)ln  - - -  (5) 

V  2  I  —  cos  V 

where  the  RMS  relative  amplitude  of  tiie  fluctuations.  For  anisotropic  fluctu¬ 
ations  the  value  of  io  in  eq.  (5)  is  replaced  by  the  value  of  4  of  the  correlation 
length  in  the  incidence  plane  in  the  direction  perpendicular  to  the  density  gra¬ 
dient,  i.e. 


9an  —  O'n 


n  I  +COsf? 

2(fc<,/x)(^of^)/n--co-^ 


(6) 


AIM  Performance 

Equations  (3-6)  were  used  to  compute  degradation  of  the  AIM  performance 
to  electron  density  fluctuations  induced  by  irregularities  in  the  ambient  neutral 
density  and  heater  control.  A  worst  and  a  best  case  were  examined.  For  the  worst 
case  self-absoiption  was  neglected  and  the  final  electron  density  was  taken  as 

n{x)  =  (7) 

where  i/(x)  is  the  local  ionization  rate  and  t  the  heater  pulse  length.  For  the  best 
case  it  was  assumed  that  self-absorption  limits  the  electron  density  so  that 

71  (z)  ~  i/"(x)  (8) 

with  a  positive  (a  «  1  —  2).  The  following  results  were  derived. 


VI 


1.  Neutral  density  irregularities 


The  upper  bound  for  the  signal  loss  due  to  irregularities  in  the  neutral  density 
derived  using  eq.  (7)  is  given 


10/  VlOOm 


CIV 


(9a) 


X  sin^O  cos^O  x 


('4 


+  cos  6 
—  cos  0 


In  eq.  (9a)  A  =  10),  the  number  of  growth  times  of  the  ionization, 

Lc(~  4  X  10^  km)  the  correlation  length  of  the  neutral  density  fluctuations  and  f 
the  radar  frequency.  For  the  saturated  case  is  given  again  by 


Cn  10  ^«^(ioOm)(lOMi/2)  (■ 


.  1+COS0\ 

X  sinrO  cos^O  x  ( In - -  ) 

V  1  -  cost? ) 

Notice  that  in  eq.  (9b)  q  si  0(1)  giving  insignificant  loss. 


(9b) 


2.  Heater  induced  irregularities 

If  the  heater  induces  density  irregularities  with  variance 

f 


lOM//; 


xln 


1  +  cos  0 
1  —  cos  t? 


(10) 


v  1 1 


where  Cq  is  the  correlation  length  of  the  heater  power  density  fluctuations. 
Equation  (10)  can  also  be  cast  in  terms  of  the  value  ap  of  the  variance  of  the 
heater  power  fluctuations.  Using  eq.  (7)  to  connect  the  ionization  rate  to  the 
density  fluctuation  that  induces  we  find  as  an  upper  limit 


‘(i^) 


10t77  7  V  \QMHz 


(Ha) 


1  +  cos  6 
1  —  cos  6 


Notice  that  ap  =  10''^  corresponds  to  — 40  dB  heater  power  density  fluctuations. 
If  we  use  eq.  (8)  which  corresponds  to  the  saturated  case  £ h  will  be  given  by 


-2.1  X 


\  1  —  cos9  J 


I  ' 

lOMHz 


(lib) 


1.  THEORETICAL  FOUNDATIONS 
1.1  General  Formulation 

We  start  with  the  wave  equation  in  a  medium  with  permittivity  6(r).  For  a 
monochromatic  wave  with  frequency  uj  it  reduces  to  the  Helmholtz  equation 

V‘^(L,i)  +  ^  0  (1) 


C‘‘ 


(2) 


For  short  wavelength  propagation,  i.e.  if  0(r)  in  the  propagation  direction  (which 
we  take  as  the  x-direction)  varies  slowly  with  respect  to  the  wavelength  l/ko.  cq. 
(1)  reduces  to  the  “parabolic”  wave  equation 

with  V’  =  *ind  Ip  a  slowly  varying  function  of  x  (Flatte,  1979;  Dashen, 

1979).  If  we  write 


^(r.)  =  «(z:)  +  ^^(l) 


(4) 


i.e.  an  average  permittivity  and  a  fluctuating  parts,  eq.  (3)  becomes 


i  di>{r) 
fco  dx 


-1 


i’iL)  -  ^[1  - 
l] 


(5) 


Equation  (5)  provides  a  clear  distinction  of  the  factors  affecting  the  range  propa¬ 
gation  of  ip{r).  The  first  term  on  the  rhs  describes  diffraction  due  to  the  transverse 


curvature  of  ip{r),  the  second  term  refraction  due  to  the  spatial  variation  of  t(r ) 
and  the  third  term  scattering  due  to  the  random  variation  Se(r).  If  we  consider 
propagation  from  a  point  source  located  at  the  origin  of  coordinates  to  a  point  r 
ipiz)  can  be  written  as 

iiL)  =  T— “(n)  (6) 

47rr 

where  the  factor  accounts  for  energy  flux  conservation  due  to  geometric 
spreading,  while  u(r)  incorporates  phase  effects  (including  absorption)  during 
propagation. 


(10) 


We  assume  next  that  the  eikonal  koSi(path)  is  a  Gaussian  variable  and 

(i)  Interchange  averages  and  integrations. 

(ii)  Use  the  fact  that  for  a  Gaussian  variable  at 

<  e'"  >=  ^  ^ 


Under  these  conditions  eq.  (10)  becomes 


<u>^  N 
xexp 


j  d{path)eip[ikoSo{p<ith)] 

(-\kl  <  SUp^th) 


(11) 


where 

<  Slipath)  >=  j  <  (  j  >  (12) 

path 

Equation  (11)  gives  a  physically  transparent  but  deceptively  simple  looking  result. 
The  first  term  describes  the  propagation  in  a  smooth  medium,  while  the  term 

R  =  exp(^-^kl  <  Siipath)  (13) 

decribes  the  modification  of  the  propagation  by  the  presence  of  fluctuations.  For 
jR  1  we  recover  propagation  through  a  smooth  medium.  From  eq.  (13)  the 
signal  reduction  C  will  be 

C  ^ -A.3[kl  <  S^iipath)  >]  dB  (14) 

In  practical  terms  the  important  task  is  the  to  determination  of  the  value  of 
<  S^ipath)  >  under  the  relevant  propagation  conditions.  In  the  next  section 


3 


we  will  consider  a  reflection  from  a  finite  size  surface.  This  will  allow  us  to 
compare  eq.  (12)  with  standard  results  (Beckman  and  Spizzichino,  1963],  The 
AIM  case  will  be  examined  in  the  subsequent  sections. 


1.2  Reflection  by  a  Finite  Size  Surface 
1.2.1  Smooth  surface 

We  apply  first  the  above  results  to  the  case  of  reflection  by  a  smooth  surface. 
Consider  a  surface  with  dimensions  ±X  and  ±Y  (Figure  1).  The  value  of  So  (path) 
will  be  given  [Ogilvy,  1991] 

So{xo,yo)  =  Axo  +  Byo-\-Ch  (15a) 

A  =  sinOi  —  sin92  cos9z 

B  =  —stn92  sin9z  (15b) 

C  =  — (cos^l  +  COs92) 

and 

<  u  >=  Ncos9i  j  dxodyo  x  explikoSo{xo,yo)]  x  Rixo,yo)  (16) 
For  a  smooth  surface  (i.e.  R  =  1)  with  area 

Am  =  4Xr  (17) 


<u>=  Ne**"^^*  x  cos9i  X  Am 

sinkoAX  sink^BY 
^  k„AX  ^  koBY 


(18) 


4 


Absorbiilg  the  constant  phase  factor  exp(likoCh)  to  the  normalization  N,  we 
recover  the  standard  formula 


<  u  >=  N(A\i  cos6\) 

sinko  AX  sinko  BY 
^  koAX  koBY 
For  specular  reflection  A  =  0.  B  =  0 


(19) 


<  u  >  —  N{Ai^i  cosd\) 

spec 


(20) 


1.2.2  Rough  surface 

Let  us  examine  next  specular  reflection  from  a  rough  surface.  From  eq.  (16) 
with  A  =  B  =  0,  we  find 


<  u  >=  N  cos6i  I  dxodyo R{xo,yo)- 

3P  J 


(21) 


For  a  surface  with  random  height  inhomogeneities  with  RMS  standard  deviation 
A/i  independent  of  the  position  (xo,  y©)  we  find,  using  eq.  (13)  in  eq.  (21), 


<  u  >=  N  cosO\  Am  eip(  — 4cos^5i 

sp  \  2 


or 


>(-|) 


<  u  >=<  u  >  eipl 

sp  sp  V  2 


go  =  4kg  A/i^  cos^O] 


(22) 


(23) 


(24) 


The  loss  during  specular  reflection  will  be 


£(rough  surface)  =  — 17(/co  A  lif  cos^O  dB 


(25) 
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The  more  general  formula  reads 


with 


<  u  >  =  <  Uo  > 


(26a) 


g  =  klAh^{cos$i  +  cos02)~ 


(26b) 


For  the  specular  case  (6i  =  ^2).  eqs.  (26)  reduce  to  eqs.  (23)  and  (24). 
1.2.3  Raleigh  criterion 


The  Raleigh  criterion  is  a  criterion  that  determines  the  degree  of  roughness 
of  a  surface  on  the  basis  of  a  simple  physical  argument.  Referring  to  Figure  2  we 
consider  scattering  in  the  azimuthal  plane  (i.e.  the  (x-z)  plane  of  incidence).  The 


difference  in  phase  between  two  scattered  rays  is 

=  ko[{zi  -  Z2){cos0i  +  COS02)  +  {^2  -  a^i)] 

x(sm0i  —  sm02) 

For  specular  scattering  (^1  =  ^2)  the  phase  difference  becomes 


(27) 


A<f>  =  2koAh  cosO\ 


(28) 


where  Ah  =  z\—  Z2.  the  interference  between  the  two  rays  depends  on  A<^.  For 
A<i)  <<  TT  the  waves  will  interfere  constructively.  For  A<i>  ~  x  the  waves  will 
interfere  destructively,  leading  to  very  weak  contribution  in  the  specular  direction. 
The  “Rayleigh  Criterion”  states  that  if  Ah  is  the  surface  RMS  deviation,  the 
surface  is  smooth  for  A</>  <  This  stated  mathematically  as 

koAh  cos0\  <  —  (29) 

4 
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From  eq.  (26b)  we  find  that 


g  =  AkzAli  cos^Oi  <  —  (30) 

4 

which  from  eq.  (25)  corresponds  to  a  loss  of  less  than  10  db,  a  rather  significant 
loss.  A  better  criterion  will  be  k^Ah  cosOi  < 


13.  Reflection  from  AIM  in  Uie  Presence  of  Fluctuations 
1.3.1  Isotropic  case 


From  eqs.(ll-13)  we  find  that  the  deterioration  of  the  signal  due  to  inhomo¬ 
geneities  in  the  AIM  surface,  will  be  given  by 
<  tx  > 


<  u  > 

O 


=  exp 


<  Slipath)  > 


(31) 


where 


(32) 


Let  us  first  examine  the  case  of  isotropic  fluctuations  with  scale  length  much 
shorter  than  the  path  through  the  cloud.  In  this  case  (Rytov  et  al.,  1990) 


<  Slipath)  >=  6al  (33) 

The  eikonal  variation  is  given  by 

^o(^)  OO 

=  f  ~~  J  Sf{k,z{s))kdk  (34) 

O  O 

where  £o(^)  is  the  zero  order  path  of  the  ray  incident  at  an  angle  0  and  Se(k,z)  the 
spectrum  of  the  fluctuations  6c  at  the  point  z.  Equation  (34)  can  be  analytically 

7 


8 


so  that 


1  +  cosO\ 
1  —  cosO ) 


The  signal  deterioration  will  be  given  by 


C  —  —4.3  g  dB 


(43) 


(44a) 


where 

^ (r-cos^) 

From  eq.  (34)  it  is  easy  to  see  that  6<t^  maximizes  when  e(z)  becomes  small,  i.e. 
near  the  turning  regions  of  the  layer.  It  is  easy  to  find  that  50%  of  the  contribution 
to  comes  from  a  fraction  of  the  layer  near  the  top  with  width  AL  given  by 

AZr  1  5m^(l  —  sin0) 

Zi  2  cos^9 


This  corresponds  to  about  a  fifth  of  the  layer  9X6  -  45°  and  1/6  at  (?  =  30°.  Notice 
that  in  this  region  the  rays  are  almost  parallel  to  x. 

A  simple  but  intuitive  way  to  derive  the  above  result  is  as  follows.  Consider 
phase  fluctuations  over  a  ray  of  length  i  and  inhomogeneities  with  scale  The 
ray  path  will  accommodate  N  ~  h  independent  inhomogeneities.  After  a  simple 
inhomogeneity  the  phase  shift  of  the  wave  wiU  be 

(jJ 

c 
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Since  the  random  quantities  corresponding  to  the  different  inhomogeneities  are 
independent,  the  total  mean  square  of  the  phase  shift  will  be 


<  52  >=  TV  <  {e,k„Sef  U]kla] 


(45) 


ikoif){koz)a'^ 


which  is  comparable  to  eq.  (43). 
13.2  Anisotropic  case 


For  a  situation  with  anisotropic  fluchiations  eq.  (34)  becomes 

^o(d)  00 

/  i(^)i  / 

o  —00 

where  refers  to  the  transverse  direction  at  every  point  on  the  path.  The  role  of 
anisotr<q)ic  fluctuations  can  be  understood  by  considering  the  value  of  the  integral 

00 

J  =  J  S,ikj_,0)d'^k^  (47) 

— C» 

for  the  case  that  the  fluctuations  are  described  by 


Ail)  =  exp 


2{a^^  cy 


(48) 


In  this  case  6e(r)  has  a  spectrum 


s,(k)  =  +  kV,  +  c^kl) 

(27r)  ^ 


(49) 


Assuming  propagation  on  the  x-z  plane  at  an  angle  7  to  the  z-axis 


ni)  = 


ac 


\/2x  yfa^cos^'y  -f  (P’sin?') 


(50) 
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Notice  that  lor  propagation  in  the  x-direction  (7  =  |), 


and  for  7=0, 


(51a) 


m  =  ^  (51b) 

Namely  the  eikonal  correlation  is  controlled  by  the  correlation  length  along  the 
propagation  direction.  From  eqs.  (46-49) 

=  (y )  X  “ 

^  r  _ds _ 1 _  (52) 

J  €(2(5)]  ^5^27(5)  +  2^00527(5) 


where  7(5)  is  the  local  propagation  angle.  For  isotropic  spectra  a  =  c  and  we 
recover  the  result  of  eq.  (41)  with  replaced  by  a.  Since  the  integral  of  eq.  (52) 


maximizes  when  c(z(s)]  is  minimum,  i.e.  near  the  turning  regions,  7(5)  w  |  and 
as  long  as  ^cos^7(s)  <<  1, 


to{e) 


o 

For  a  linear  density  profile  eq.  (53)  gives 


ds 


1^(5)] 


(53) 


(54) 


The  result  is  similar  to  eq.  (42)  but  with  the  isotropic  scale  length  £„,  replaced 
by  the  scale  length  a  in  the  x  direction.  The  value  of  g  is  given  by 


cosO 


(55) 
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2.  SOURCES  AND  SPECTRA  OF  PERMITTIVITY  FLUCTUATIONS  IN  AIM 


For  practical  AIM  systems  there  are  two  sources  of  permittivity  fluctuations 
Se(r).  due  to  electron  density  fluctuations  within  the  cloud.  Both  are  associated 
with  fluctuations  in  the  local  ionization  rate.  The  first  source  is  due  to  natural 
fluctuations  in  the  ambient  neutral  density.  The  second  is  due  to  spatial  variation 
of  the  incident  RF  power  caused  by  scanning  and  heat  control  errors.  In  both  cases 
the  electron  density  fluctuations  are  induced  by  fluctuations  in  the  ionization  rate 
i/(x).  In  deriving  the  spectrum  Sn(k)  from  the  spectrum  of  the  ionization  rate 
fluctuations  Si/(k),  we  use  two  models.  An  upper  limit  on  the  value  of  Sn(k) 
can  be  found  by  ignoring  self-absorption  in  the  formation  of  the  AIM  cloud  and 
assuming  that 

n(i)  =  (56) 

where  r  is  the  pulse  length  of  heater  and  n^,  the  initial  electron  density.  A  more 
realistic  estimate  can  be  found  by  considering  the  saturated  case  (i.e.  including 
self-absorption)  in  which  case 


n(x)  ~ 


(57) 


with  a  «  1  —  2.  The  spectrum  Sn(k)  is  defined  here  in  a  normalized  form,  i.e. 

_  f  <nn  I  I  > 

"  (2^)^  J 

so  that 


(58) 


Snik)  =  S,{k). 


(59) 
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2.1  Unsaturated  Case 


For  the  unsaturated  case,  we  assume  that  the  density  at  a  point  x,  is  related 
to  the  ionization  rate  by  eq.  (56).  We,  furthermore,  normalize  i/(x)  to  the  pulse 
length  T,  so  that  i/(x)  =  t/(x)T.  As  a  result 

<  n(i)n(i')  >=  nl  <  >  (60) 

Writing 

i/{x)  =  17{x)  +  8i'{x)  (61) 

we  find 

<  n(i)n(i')  >=  >  (62) 

Assuming  that 

<  Si/{x)  >=  0  (63) 


with 

\  <  > 

we  find  that  eq.  (62)  becomes 


<  >=  eip 


(64) 


<  n(a:)n(a:')  >=  ^exp[<  Sv{x)8u{x)  >]  (65) 

As  a  result 

<  n(x)n(x')  >  —  <  n^(x)  >  =  exp(<  iP'  >  ) 

(66) 

x[exp(<  6l'{x)8i/{x')  >)  —  1] 
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If  we  define 


/\  _  <  n(i)nv.T  )  >  -  <  n^(i)  > 

G(i  -  I  1  =  - y— - 

^  ’  <  rO-{x)  >  (67) 

=  (eip(<  >)  —  1] 

the  one  dimensional  electron  density  fluctuation  spectrum,  will  be  given  by 


(68) 


or 

Snik)  =  Jdx  -  1 

Returning  to  dimensional  vaiiables  and  defining 


(69) 


A  =  1/7 


(70) 


and 


n  ,  I  ^  >  2  /  ^ 

5„(x)  =  -  - - =  <7^F^{x) 

V 


(71) 


with  F„(i  =  0)  =  1,  F„(x  — >  oo)  =  0,  eq.  (69)  becomes 


Sn{k)  =  ^  e-‘"{exp[A2<7;F,(i)]  -  1}  (72) 

Notice  that  the  value  of  is  the  strength  parameter,  since  F„(i)  <  1.  We 
distinguish  the  cases  of  weak,  X^al  <<  1,  and  strong,  A^cr^  >  1,  fluctuations. 
2.1.1  Weak  fluctuations 


We  define  as  weak  fluctuations  if  «  1.  For  example,  if  A  ^  10  weak 
fluctuations  Imply  ^  <  .1.  The  practical  applications  require  such  as  a  condition. 
From  eq.  (72)  we  find 

Sn{k)  =  X'^S^ik)  (73) 
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Namely  the  spectral  distribution  of  the  density  fluctuations  will  be  similar  to  the 
ionization  rate  fluctuations,  but  their  strength  will  increase  by  A^. 

2.1.2  Strong  fluctuations 

In  this  case  X^al  >>  1  and 

Snik)  =  ^Jdx  e-*'--"  exp{\^alF,{x))  (74) 

In  order  to  proceed  we  assume  that 

f.(x)  =  pIjj  (75) 

This  is  a  Lorentzian  spectrum  with  correlation  length  b.  For  F|/(x)  given  by  eq. 
(75)  the  integral  in  eq.  (72)  can  be  performed,  using  the  technique  of  stationary 
phase.  We  find 

/tt  1  (  \ 

The  importance  of  this  result  lies  in  the  fact  that  while  for  weak  fluctuation  there 
is  essentially  no  spectral  shift,  for  strong  fluctuations  there  is  a  complete  shift  of 
the  spectrum  towards  shorter  wavelengths.  The  equivalent  correlation  length  is 
smaller  by  ^  and  can  be  strongly  reduced  by  adjusting  A  to  large  values. 
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2.2  Saturated  Case 


For  the  saturated  case  the  density  is  connected  to  the  ionization  rate  via  eq. 
(57).  In  this  case 

<  n(x)n(i')  >  =  <  > 


=  <  (i/(x)  +  6i/(i))“(i/(x')  +  (5i/(x'))“  > 
(-(  112^*  A  ,  2<f>v{x')Sv{x)> 

=  (i/(i))  1+0  — 


From  eq.  (77)  it  is  straightforward  to  find  that 


(77) 


Sn{k)  =  a'^S,{k) 


(78) 


Namely  the  spectrum  of  the  fluctuations  remains  the  same  but  the  amplitude  is 
modified  by  the  factor  o^. 


3.  DEGRADATION  OF  AIM  PERFORMANCE  DUE  TO  FLUCTUATIONS 
IN  THE  LOCAL  IONIZATION  RATE 

3.1  Fluctuations  of  the  Ambient  Neutrals 

The  ionization  rate  is  a  linear  function  of  the  neutral  density  N.  As  a  result, 
if  all  other  factors  are  constant 


Su{k)  =  Sjv(k)  (79) 

Based  on  available  measurements  [  Jursa,  1985]  for  the  relevant  altitude  range  of 
60-80  km,  SN(k)  can  be  described  by  an  isotropic  spectrum 

SnU  ^  I)  =  ^  (80) 
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This  corresponds  to  an  isotropic  spatial  autocorrelation  function  <  6NS.\'  |  r  > 
which  for  “small”  distances  has  a  linear  decay  with  radius  with  scale  length 
Lc  —  The  wavenumber  spectrum  of  eq.  (79)  has  a  spectral  index  of  p=4, 
i.e.  decays  as  The  corresponding  one  dimensional  wavenumber  spectrum 
has  a  spectral  index  o  =  p-2  =  2.  From  Figure  3  and  Table  1  we  see  that 
Lc  ~  4  X  lO^kni,  giving  kc  ~ 

For  the  AIM  problem  the  maximum  size  of  the  inhomogeneities  of  interest 
corresponds  to  the  minimum  of  either  the  Fresnel  size  Lp,  which  is  of  the  order  of 
a  few  kms,  or  of  A(<?)  defined  by  eq.  (39).  It  is  convenient  to  cast  eq.  (79)  as 

SN{k)  =  ^  (81) 


where 


A  = 


(82a) 


k 


m 


(82b) 


From  eqs.  (73),  and  (82) 


s,{k)  =  S„{k)  =  k  >  kr„ 


(83) 
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For  the  linear  profile  eqs.  (40)  and  (82)  give 


OO 


(84) 


For  the  general  AIM  case  A  <  Lp  and  As  a  result 


From  eqs.  (14)  and  (34)  the  signal  loss  will  be  given  by 


(85) 


C  =  -A.dirX^{koL){k„  A  (6)) 

,  1  +  cose 

xln - -  dB 

I  —  cose 


m 


(86) 


For  a  linear  profile  using  the  results  of  Section  1 .3,  eq.  (86)  becomes 
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'  ■  - '(a'(ii9‘(B75;)’{^)  -  ('ra) 

(87) 


xsin^O  cos^O  (IB 

Equation  (87)  represents  an  upper  limit  in  the  fluctuations,  since  as  it  can  be  easily 
verified  in  this  case  the  weak  fluctuation  limit  is  always  valid.  For  the  saturated 
case  the  loss  will  have  the  same  form  as  in  (86)  but  with  A  replaced  by  a.  Since 
a  0(1)  the  loss  rate  will  be  totally  negligible,  in  this  case. 


3.2  Effect  of  Heater  Induced  Fluctuations 
3.2.1  Loss  ^  a  function  of  Um  density  fluctuation  level 

Fluctuations  in  the  electron  density  can  also  arise  from  the  finite  step  of  the 
heater  during  painting  of  the  AIM  cloud,  as  well  as  from  any  other  heater  phase 
errors  [Short  et  al.,  1990].  The  analysis  below  addresses  the  constraints  imposed 
by  a  heater  step  size.  The  analysis  relevant  to  this  case  is  the  anisotropic  analysis 
discussed  in  Section  1.3.2.  Assuming  that  the  heater  induced  fluctuations  have  a 
correlation  length  a  in  the  ray  direction  equal  to  the  step  size  slq  of  the  ionizing 
beam,  we  find  that  the  power  loss  will  be  given  from  eqs.  (25)  and  (56)  as 


For  a  step  size  of  10m,  f  =  10  MHz  and  L  «  100m,  the  round  trip  loss  will  be 
4-5  db  for  10%  density  fluctuations.  To  achieve  the  same  performance  with  50 
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MHz  the  fluctuation  level  should  be  reduced  to  2%.  As  discussed  by  Short  et  al. 
[1990]  these  can  be  easily  achieved 

3.2.2  Loss  as  a  function  of  heater  power  fluctuations 


In  this  case  the  fluctuations  will  be  anisotropic  and  we  should  use  the  equations 
derived  in  Section  1.3.2  to  calculate  the  loss.  We  assume  that  the  power  fluctuation 
spectrum  is  given  by 


S,(i)  =  (89) 

and  taking  i/  ~  P  we  find 

5„(x)  =  (90) 


For  the  case  of  weak  fluctuations 


«  1 


(91) 


5„(x)  =  AVjFp(^) 

For  a  Gaussian  spectrum  Fp(|),  eq.  (56)  in  conjunction  with  eq.  (12b)  gives 

1  I _ /I  ^  ' 


\+cos0 

xln - — -  dB 

1  +  Sind 


Notice  that  ap  =  10"^  corresponds  to  power  fluctuations  of  4-40  dB.  According 
to  Short  et  al.  [1990]  power  fluctuations  as  low  as  — 60  dB  can  be  achieved. 
For  the  saturated  case  replacing  A  by  o  we  find  losses  by  at  least  an  order  of 
magnitude  smaller. 


20 


3.2.3  Strong  fluctuations 


It  is  interesting  to  examine  the  upper  limit  of  loss  for  the  case  of  strong 
fluctuations  in  the  sense  of  Acrp  »  1.  Such  is  the  case  if  the  heater  power 
fluctuations  are  large  (i.e.  >  10  db).  In  this  case  using  eq.  (46)  for  Sa^,  with  the 
spectrum  given  in  eq.  (77)  we  find 

(r^) 

Notice  that  the  correlation  length  has  shifted  to  instead  of  b.  The  loss 

Up  A 

computed  on  the  basis  of  eq.  (91)  is 
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Figure  1. 
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Figure  2. 
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Table  1.  Estimated  rms  Differences  (%  of  mean)  Between  Densities  at  Locations  90, 
1 80,  and  360  km  Apart  During  the  Midseason  Months  in  the  Tropics 


Alt 

(km) 

90 

minm 

133 

360 

90 

6)1 

360 

90 

wmm 

■Tfil 

360 

90 

Oct 

180 

km 

360 

10 

0.10 

0.13 

0.18 

0.10 

0.13 

0.18 

0.10 

0.13 

0.18 

0.10 

0.13 

0.18 

15 

0.13 

0.17 

0.25 

0.11 

0.14 

0.21 

0.16 

0.20 

0.30 

0.16 

0.20 

0.30 

18 

0.50 

0.61 

1.00 

0.34 

0.42 

0.68 

0.30 

0.37 

0.60 

0.34 

0.42 

0.68 

20 

0.28 

0.34 

0.56 

0.28 

0.34 

0.56 

0.24 

0.29 

0.48 

0.24 

0.29 

0.48 

25 

0.28 

0.34 

0.56 

0.28 

0.34 

0.56 

0.24 

0.29 

0.48 

0.26 

0.32 

0.52 

30 

0.30 

0.37 

0.60 

0.30 

0.37 

0.60 

0.28 

0.34 

0.56 

0.30 

0.37 

0.60 

35 

0.34 

0.42 

0.68 

0.30 

0.37 

0.60 

0.30 

0.37 

0.60 

0.36 

0.44 

0.72 

40 

0.40 

0.49 

0.80 

0.44 

0.54 

0.88 

0.48 

0.59 

0.96 

0.44 

0.54 

0.88 

45 

0.46 

0.56 

0.92 

0.40 

0.49 

0.80 

0.60 

0.73 

1.20 

0.52 

0.64 

1.04 

50 

0.56 

0.69 

1.12 

0.54 

0.66 

1.08 

0.72 

0.88 

1.44 

0.54 

0.66 

1.08 

55 

0.66 

0.81 

1.32 

0.56 

0.69 

1.12 

0.84 

1.03 

1.68 

0.78 

0.96 

1.56 

60 

0.84 

1.03 

1.68 

0.66 

0.81 

1.32 

1.00 

1.22 

2.00 

0.82 

1.00 

1.64 
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